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1. Introduction 

A quantum system with continuous degrees of freedom can be represented in terms of a 
Wigner Function [2] denned as a real function on the phase space. The Wigner function 
(WF) is similar to a probability distribution (its integration over the phase space is 
normalized to one) even if it can take negative values on restricted domains. There is 
an extensive literature on the continuous- system WF (31 -0 due to its wide applicability 
in different contexts of physics. Concerning the WF of discrete systems, the literature 
is less extensive even if this theme has recently attracted a lot of interest mainly in view 
of the role that a discrete phase-space structure can play within Quantum Information 
Theory [JJ, [H]-[II|- In this respect paper j7] contains a useful list of references. Different 
generalizations of the WF to quantum systems with a finite-dimensional Hilbert space 
have been proposed in the literature such as 1) the continuous WF for spin variables 
[To] and 2) the definitions of WF based on a discrete phase space. As to the latter, early 
studies were made in [TOl UZ| . A discrete WF has been introduced in (Tj and |TH] which 
generalizes the 2x2 case of J7j and is valid for systems having a iV-dimensional Hilbert 
space, with N a prime or a power of a prime. More recently, an alternative definition 
of WF involving Galois fields [H|-[B] has allowed the study of the composite-dimensional 
case and evidenced several interesting tomographic properties [HI HI] • 

The present article is focused on studying the entanglement properties of a two- 
qubit system by using the discrete WF defined in [TJ . It is worth noting that two-qubit 
(and more in general many-qubit) systems have received an increasing attention, not 
only within spin models, but also in the recent literature on optically trapped bosons 
modelled within the Bose-Hubbard picture [2D]- The impressive experimental progress 
in controlling the spatial trapping of bosons makes the realization of many-qubit systems 
a quite realistic objective. In this paper we find some new separability criteria and recast 
other known criteria in terms of discrete WF's. We check them evaluating two-qubit 
entanglement and show that the discrete WF describes both classical and quantum 
correlations better than the density-matrix approach. We note how using the WF not 
only improves the visualization of the system state but is also expedient experimentally: 
since the WF is directly related to tomographic techniques, the separability criteria 
coming from the WF do not require to know of all the matrix elements. 

In section 2, we review the definition of discrete phase-space given in reference [TJ 
-this is particularly useful for our purposes- and the corresponding discrete Wigner and 
Characteristic functions. We present some basic properties of these functions which 
provide a useful tool for studying quantum correlations. In section 3 we consider 
four different separability criteria in terms of the two-qubit discrete WF and of its 
covariance matrix. In particular, we show that 1) there is a negative value of discrete 
WF that allows one to discriminate between separable and entangled states, 2) the 
Partial Transposition criterion can be reformulated in terms of the two-qubit discrete 
WF, 3) there is a nontrivial link between the separability of the density matrix and 
the non-negativity of the WF's corresponding both to the density matrix and to its 
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partial transposed matrix, 4) the Local Uncertainty Relations relevant to phase space 
operators can be generalized in terms of the WF covariance matrix (thus evidencing 
the difference between classical and quantum correlations), and 4) the Generalized 
Uncertainty Principle, so far studied for WF's relevant to continuous phase space, is 
extended to the case of a discrete WF. 

2. Discrete phase space and Wigner function 

In a discrete r-dimensional Hilbert space, with r a prime number, the phase space can 
be defined pQ as a r x r array of points. The latter can be labelled by pairs of coordinates 
a = (q,p), each taking values from to r — 1. For each coordinate we define the usual 
addition and multiplication mod r thus obtaining the structure of a finite mathematical 
field F r with r elements (0, 1, r — 1). If the dimension is N = r n , with r prime and n 
an integer greater than 1, the discrete phase space can be built in two ways, both giving 
a discrete phase space formed by a, N x N grid: the first involves the extension F^r of the 
primitive field F r [ZJ E] , while the second is based on performing the n-fold cartesian 
product of r x r phase spaces pQ . In the present article we will use this last definition of 
discrete phase-space (entailing the definition of WF given in pQ). The choice of phase- 
space structure is justified by the direct connection with the tensor-product structure 
of the Hilbert space ensuing from the decomposition of the system in two or more 
subsystems, which is a useful feature for studying the entanglement. According to this 
definition, the phase-space points a are labelled as n-tuple (ati, a 2 , a n ) of coordinates, 
each ai pertaining to the z-th subsystem. In each subsystem with prime dimension r we 
can build standard lines as set of points satisfying equation (uq + vp) mo dr = c. However, 
we cannot define uniquely lines over the entire phase space (with modular arithmetic): 
in reference |Zj there is an example of two sets of points which form two parallel "lines" 
but intersect in two distinct points. Nevertheless, we can define the alternative concept 
of slice PQ: given a set of n lines {A,} (one for each subsystem), the slice is the set of 
all points a = («i, a 2 , a n ), where G Aj. A weaker notion of parallelism can be 
defined: two slices are parallel if each of the n lines forming the first slice are parallel 
to the corresponding n lines of the second slice. 

2.1. Definition of the discrete WF 

The discrete WF relevant to a iV = r"-dimensional system can be defined PP by means 
of the set of discrete phase-point operators A(a) (or A(a) ). Consistently with 
the definition of phase space in terms (of cartesian product) of constituent subspaces 
PQ, phase-point operators A(a) are defined as tensor product of phase-point operators 
relevant to the corresponding subsystems: A(a) = A(ai) ® A(a 2 ) <8> ... <8> A(a n ). Since 
they form a complete orthogonal basis for the Hermitian N x N matrices, any density 
matrix can be written as p = J2 a W(a)A(a), where the real- valued coefficients 



W(a) 



1 



tr[pA(a)} 



(1) 



N 
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represent the discrete Wigner function (also called the discrete Weyl symbol). Phase- 
point operators exhibit two basic properties: i) for any couple of points (ai,a 2 ) 

tr[A(ai)A(a 2 )} = N5(ax, a 2 ) , (2) 

ii) given any slice A in the phase space, the projector relevant to A can be written as 
^ 1 



N 



5>( 



Oil 



(3) 



The latter definition implies that the set of all Py for which A' is parallel to A 
forms a set of mutually orthogonal projection operators. Moreover, given the slice 
A = (Ai, A2, A n ), then P\ is the tensor product P\ l (g)P\ 2 ®..®P\ n of projectors relevant 
to the subsystems. Such properties, analogous to those characterizing continuous phase- 
point operators [T], can be used to derive the discrete- WF properties. Owing to formulas 
El and El discrete WF's feature two crucial properties. First, if W(a), W(a) correspond 
to density matrices p, p', respectively, then formula 121 entails that 

Nj2W(a)W'(a) = tr(pp') . (4) 

a 

Second, due to equation El given a complete set of N parallel slices, for each slice A, 
the N real numbers p\ = J2 a ^x W a are the probabilities of the outcomes of a specific 
measurement associated with A. Hence J2 a W a = 1 (normalization property). 

Let us consider first the simple case N = 2 (single qubit). The phase-point operators 
can be written in terms of Pauli matrices as jl] 

1 
2 



A(a) 



I+(-lYa z + (-l) p a x + 



-1 



\q+p 





' 


1 " 




' -i " 


, o- x = 


_ 1 


_ 




i 



and o z = [a x , a y ]/2i. The single-qubit phase space is the set of points a = (q,p), where 
q,p = 0,1, exhibiting properties of F 2 . Since the density matrix for a general one- 
qubit state can be written (within the standard computational basis) in terms of three 
independent real elements p 00 , Re(p i), im(p 01 ), we have W(q,p) = 4{[1 + (— l) 9 ]p o + 
[1 - (-l) 9 ]pn + (-l) p 2i?e(p i) + (-l) q+p 2Im(p 01 )}. In the case of two qubits, the WF 
has a more complex expression. Upon noting that the phase space operators are defined 
as A(a±, cx.2) = A(cti) ® A(a 2 ), the WF becomes 



W(q 1 ,q 2 ,Pi,p 2 ) 



-tr\pA(qi,pi) <S)A(q 2 ,p 2 )] 



(5) 



When necessary, we shall write W p , where the subscript means that the WF is associated 
to density matrix p. 



2.2. The discrete characteristic function 

The set {J, o~ x , a y , a z } forms an orthogonal basis for the set of hermitian operators 
acting on a single qubit. Thus any density matrix p for a single qubit can be written as 
P — I v)S(u, v) while the characteristic function for a single qubit is 

X(u,v) = tr[pS(u,v)] (6) 
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where S(0, 0) = I, 5(1,0) = a x , S(0, 1) = o~ z , and S(l, 1) = oy When necessary, we 
write the argument /3 instead of (u, v), or the single index i, where i = u + 2v assuming 
integer values from to 3. In the case of two qubits, any density matrix can be written 
as P — i S/3i/3 2 x(Pii P2)S(Pi) <8> S(02)- The two-qubit characteristic function is thus 
defined as xiPiife) = t r [pS(Pi) ® ^ifo)]- Function x(P) is connected with the discrete 
WF by a discrete Fourier transform. For example, in the single-qubit case (r = 2), 
W(q,p) = \ Y J p{'~^Y qu+pvS> x{ u ^ v ) ) an d in the two-qubit case 

W{q 1 ,q 2 ,p 1 ,p 2 ) = ± X! (-l) (9lUl+Pl,,l)+(92M2+P2l ' 2) x(wi,W2^i^ 2 ) • 

It is worth noting that the determination of x{P) is connected to a specific tomographic 
technique [0]. For example, in the case of spin | particles, it consists of repeated 
measures of spin "up vs down" along three directions relevant to each particle. The 
WF can be determined via the previous equations from the characteristic function. The 
Inner Product Rule for x(P) corresponding to formula |U is 

^ExW(P)=tr(PP')- (7) 

Operators S can be thought as translation operators. The single-qubit WF relevant to 
p' = S(a, b)pS(a, b)^ is W p '(q,p) = W p (q + a, p + b) , where the sum in the argument is 
mod2. In the two-qubit case, operators S act as translation operators on each single- 
qubit phase space. This formalism supplies a useful tool for recognizing the translational 
covariance of the WF PITU]. 

2.3. Graphical representation, pure and mixed states. 

The phase space of the single qubit is represented in tableHJ left panel. A simple example 
of WF is given by the spin j = 1/2 coherent state |£) = (|0) + £ |1)) /yl + |£| 2 (where 
|0) and |1) are the spin-up and the spin-down state, respectively), which is known to 
be the most general pure state for a single qubit. Given the corresponding WF is 
depicted in tabled right panel. The WF can have only one negative element (otherwise, 
at least one probability associated to a direction would be negative) determined by 
\Re(£) +Im(£)| > 1 or \Re(£) — ira(£)| > |£| 2 . An important observation is that, while 
most positive value of the single-qubit WF is 1/2, the most negative value (1 — v / 3)/4 
is assumed in correspondence to state |£) with £ = (1 + i)/(l — [H]. These extremal 
values are useful to write a simple separability criterion. 

In table EJ we illustrate W a = W(qi, q^,Pi,P2) on the discrete phase-space points 
for the two qubit case, where the phase-space label is a = (oti, a.2) = (qi, Q2,Pi,P2)- This 
table is useful to clarify the notation we adopt for the WF (which differs from that of 
reference j^j). The purity character of a state can be evinced both from the WF and 
from the characteristic functions. From the general equations |U and we find 

1 < N £ W{af = i X(P) 2 = Hp 2 ) < 1 , (8) 
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Table 1. Left panel: graphical representation of the discrete WF for one qubit. Right 
panel: an example of WF for SU(2) coherent state with j — 1/2 (up to a factor 

1/2(1 + iei 2 )) 



1 

p 


W(0,1) W(l,l) 
W(0,0) W(1,0) 




1 

q 



1 - i?e(0 - Jm(0 


|e| 2 -i?e(0 + /m(0 


1 + Re(0 + Jm(£) 


|e| 2 + Re(0 ~ lm(0 



Table 2. Graphical representation of the discrete WF for two qubits 
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10 


W(00,10) 


W(01,10) 


W(10,10) 


W(ll,10) 


(Pl,P2) 01 


W(00,01) 


W(01,01) 


W(10,01) 


W(ll,01) 


00 


W(00,00) 


W(01,00) 


W(10,00) 


W(11,00) 




00 


01 


10 


11 






(91,92) 







where the equality holds for pure states whereas the inequality is involved by mixed 
states. In order to define a "mixed state" we recall that, upon introducing the basis 
{|e) : p\e) = p e \e)} relevant to a given density matrix p, a state is "mixed" when 
more than one eigenvalue p e is nonzero. In this case the system state is represented 
by p = S e Pe|e)(e| where probabilities p e evidence the characteristic lack of information 
about the relative phases of the state superposition. An interesting feature of formula 
IH1 is that the pseudo-probability WF can not be concentrated in a too small region of 
phase space. We will see that this is equivalent to the Uncertainty Principle. 



2.4- Axis operators 

In the single-qubit case the axis operators are defined as [T] & = \ J2qp£i(QiP) A(q,p), 
where £i(<?,p) := p, &(q,p) ■= (q + p)mod2, and &(q,p) := q while «fi = p, «f 2 = d, and 
^3 = q, (relevant to vertical, diagonal and horizontal lines) are the vertical, diagonal and 
horizontal axis operators, respectively. The explicit form of operators £j reads 

i = \[I-S(i)]. (9) 

where S(i) is defined after formula El In this case operators q and p play the role 
of (discrete) position and momentum operators, respectively. The spectrum of such 
axis operators is completely determined by the (two-eigenvalue) spectrum of a z and 
a x , respectively. In the sequel operator d will be named diagonal- direction (or simply 
diagonal) operator, since it is connected to the diagonal lines. Notice that obey 
commutators = 2ieijk£k showing the SU(2) algebraic structure. Thus they have 

essentially the same physical meaning of the three Pauli matrices which is to describe 
two-level systems. The only difference is that the relevant eigenvalues are and 1 (rather 
than ±1/2) that are more useful for treating quantum information applications. 
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Analogously to the continuous case, we can define the anticommutator 

{§,§}* = i(S6+SS) (io) 

(where the label S stands for standard). Differently from the continuous case, the mean 
value tr(p{l;i, £j}g) cannot be written as sum over the phase-space points of qpW(q,p). 
We thus introduce an alternative definition of anticommutator 

{§, §}d = 2 (S + & - l e ufclffc)> (11) 
where D stands for discrete. This definition allows one to express the symmetrized 
product {quPj}D as a sum over the phase space of W(q,p) multiplied by qp. In general, 

(£) = £6W( ? ,p), ({L£ 3 } d ) = J2^ w ^p)- (i 2 ) 

9,P 9,P 

The introduction of the symmetrized product [TT] is motivated by the identity = 
§[& + - (& + tj)modi] ( with & = q,p,de 0, 1). 

In the two-qubit case (iV = 4) we can perform nine possible measurements (nine 
combinations of Pauli matrices), corresponding to the nine striations of phase space. 
This tomographic scheme is not the most efficient since five orthogonal measurements 
suffice to determine the state. Nevertheless, we will consider such scheme (involving 
nine striation operators £j <S> £j) in that it leads to a definition of the WF exhibiting 
more interesting entanglement properties. 



3. Entanglement properties in two qubit systems 

Given a two-qubit density matrix p, such a state is said to be separable if there exists 
a decomposition p = J2kPk iV'fc) (V'fcl ® \4>k) (4>k\ (with the probabilities J2kPk = 1)- A 
nonseparable state is said to be entangled. If a state can be written as a density- 
matrix product p = p' ® p" (p', p" relevant to the two constituent subsystems) then the 
corresponding WF is written as W(qi,q2,pi,P2) = Wpt(qi,px)Wpii (92,^2)- Thus the WF 
associated to a separable state is 

W{q 1 ,q 2 ,p 1 ,p 2 ) = Y,PkW' k {q 1 ,p 1 )Wl(q 2 ,p 2 ) . (13) 

k 

In this perspective -so far scarcely considered in the literature- an entangled state 
exhibits a WF that cannot be written in the form ED If a classical probability 
distribution can be written non-trivially as p(qi,q 2 ,Pi,p 2 ) = J2iPkPk(QiyPi)Pk( ( l2,P2), 
the presence of more than one pk > indicates a (classical) correlation. The 
WF representation [T3I clearly evidences that separable states display a classical-like 
correlation since the related WF's have the same form of a classical distribution, whereas 
entangled states embody a different type of correlation named quantum correlation. 
We investigate the entanglement properties of two-qubit WF within 1) the negativity 
approach 2) a direct reformulation of PT criterion in terms of WF, 3) the study of 
non-negativity of WF relavant both to the density matrix and to its partial trasposed 
(deriving from the PT criterion) 3) the Local Uncertainty Relation (LUR) approach and 
4) the Generalized Uncertainty Principle (GUP) of the continuous case. 
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Table 3. Left panel: graphical representation of the discrete WF for a two-qubits 
separable state with the most negative values. Right panel: graphical representation 
of the discrete WF for the singlet state 
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singlet WF 



3.1. Negativity of WF and entanglement 

We show that the negativity of W p can be connected to the non-separability. We give 
a sufficient condition for non-separability, based on the observations of subsection 12.31 
where it is shown that any single-qubit WF assumes (1 — Vo)/8 as most negative 
value and 1/2 as most positive value. We can get a two-qubit WF with negative 
elements considering the product of WF's of single qubit W \ai)W (a 2 ) , where W{a\) 
has negative elements, while W(a 2 ) is positive. The most negative value of such a 
two-qubit WF is given by considering the most negative value for iy(ai) and the most 
positive value for W(a 2 ), as exemplified in table El left panel. The minimum value 
we get is (1 — Vo)/8 ~ —0.0915, which is the lower limit not only for WF relevant to 
product states. It is easy to show that any convex combination [T3] (i.e. separable states) 
have, as most negative value, (1 — \/3)/8. However, the value we have found is not in 
general the most negative value of a WF, as we can see in table H3 right panel. The 
state represented is the singlet state (a particular the Bell state), which results to be 
maximally entangled. Thus if the WF has a negative value H^(a) < (1 — \/3)/8, the 
state is entangled. Of course, it a WF has all the values W(a) > (1 — v / 3)/8, then the 
state can be entangled or separable. The Partial Transposition criterion, analyzed in 
next two sections, will be useful in such cases. 



3.2. Partial transposition criterion 

For all bipartite states (both discrete and continuous), the well-known partial- 
transposition (PT) criterion (23 EOj turns out to be a necessary condition for 
separability. In the 2x2 and 2x3 dimensional cases, it is also a sufficient condition. 

In the discrete case, the transposition action on a single-qubit WF and on its 
characteristic function gives respectively W p r(q,p) = W p (q,p) — (—l) q+p tr(pa y ) and 
X p t(1, 1) = — x p (l, 1) (where x is unchanged for (u,v) ^ (1, 1)). In the two-qubit case, 
the PT with respect to the second subsystem of p provides the new operator p T2 whose 
matrix elements are Pmimu = Pmunfj,, where latin (greek) indices refer to the first (second) 
subsystem. The WF W p T 2 (a.\, a 2 ) corresponding to p T2 reads 

j Pmunfi A nm (ai) A^(a 2 ) = | ^ p mi mvA mn (oti) A*^ v {a 2 ) = jtr[pA(ai) ®A*(a 2 )] , 
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Table 4. Application of the PT criterion to a Werner state. 
W p (a) W p T(a) 
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where partial transposition is shown to be equivalent to the substitution A ai <S> A a2 — > 
A ai (g> A* 2 . Interestingly, the latter is connected to the alternative definition [HJ d of 
WF involving tomographic properties that differ from those discussed in section E31 As 
to the action of PT on a WF and its characteristic function we find 

W pT2 (a) = W p (a) - T p (a) , x P ^ H) = -Xp(«i, n ) > ( 14 ) 

where the trace-like term T p (a 1 ,a>2) = (—l) q2+P2 tr(pA ai g) a y )/2 in equations flU 
embodies the effect of the PT. It is known that the operator p T2 relevant to a separable- 
state density operator possesses non-negative eigenvalues. In view of the properties just 
discussed, the PT criterion can be reformulated within the WF approach. If p T2 has all 
nonnegative eigenvalues, then tr(p T2 p') > 0, for all density matrices p', thus giving 

Y,WM a ) W A a ) =E^(«)Xp'(a) >0,W^(a),v(a), (15) 

a a 

which is a necessary and sufficient condition for separability. To illustrate this result, 
we consider the Werner (mixed) state p = x\ty~)(^~\ + (1 — x)I/4, where \^~) = 
(|0, 1) - (1,0))/^ and assume p' to be the pure state |$+) = (|0,0) + |1, 1>) /a/2- It is 
easy to show that (see table EJ) J2 a W p r 2 (a)W$+(a) = (1 — 3x)/16, consistent with the 
well known separability of the Werner state for x < 1/3. This method states that the 
separability of p is ensured when inequality E3 holds for any W' a thus having a limited 
operational value. Nevertheless, it is important in that 1) its violation for some W' a 
entails that p is entangled, and 2) it is useful to link in a direct way the non-negativity 
of the WF to the entanglement properties, as we show in the next section. 



3.3. Non-negativity of WF and separability 

We now return to the difficult problem of establishing a connection between the 
non-negativity (negativity) of WF and the separability (non-separability) of the 
corresponding state. The starting point consists in observing that any Bell state has 
a WF with negative elements, whereas a Werner state has non-negative WF for any 
x < 1/3 (separable cases) (this is illustrated in table EJ). In this respect, however, 
we know that exist separable states with negative WF such as the state in table 13.11 
left panel. On the other hand, one might conjecture that the non-negativity of WF is 
a sufficient condition for separability. Unfortunately, one can show that non-negative 
WF's exist which correspond to entangled states. As a possible strategy for solving 
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Table 5. Two examples of covariance matrix of two qubit WF 
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this problem, we thus propose a simple method based on considering the non-negativity 
features of the WF's relevant both to p and to p T ' 2 to check the separability of the state. 

Let us assume that, given a state p with W p and W p r 2 no n- negative for every phase- 
space point a, there is a W p i giving J2 a W p T 2 (a)W p '(a) < (this entails, using H31 that 
the state is entangled). We show that these assumptions lead to a contradiction. First 
we note that the last inequality can be rewritten as J2 a W p (a)W p '(a) < J2 a T p (a)Wp?(a). 
On the other hand, the non-negativity of both W p and W p r 2 , and equation O [which 
gives W p (a) > r p (a) for all a] imply that J2a W p (a)W p >(a) > J2 a T p (a)W p >(a), which 
clearly involves a contradiction. 

It follows that, given a state p, if both W p and W p r 2 have non-negative elements, 
than p is separable. Viceversa, if a state p is entangled, then W p or W p r 2 has negative 
values. Such a result -which is a necessary condition to ensure entanglement (sufficient 
condition for separability)- relates the nonclassic character of entangled states to the 
presence of negative elements in W p and W p r 2 . This criterion has been confirmed by 
testing it on thousands of randomly-generated density matrices and on the Werner state 
(for x < 1/3 both W p and W p r are positive, which implies separability). 

3.4- Local Uncertainty Relation 

It is known that the violation of local uncertainty relations (LUR's) is a signature of 
entanglement [23 EE]- Given two qubits 1 and 2, the inequalities 

E^[S 1} ]>^ ^d ££/[£P]>~, (16) 

i i 

are known to be uncertainty relations relevant to the single qubit systems k = 1,2, where 
= ((fi fe ' > ) 2 ) ~ (ii^') are the uncertainties relevant to the set of axis operators i} k ^ 
defined in formula 01 We note that simple calculations prove the equivalence between 
formula El and El entailing that the pseudo-probability can not be concentrated in a 
too small region of the phase space. As shown in j2Zj, in the two-qubit case, separable 
states are constrained by the single-qubits uncertainty relations 

E» r [2 1) ®J+^®8 a) ]>i- ( 17 ) 

i 

As a consequence, a state appears to be entangled if inequality H3 is violated. 
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LUR inequalities can be formulated in terms of WF, by defining the first-order 
covariance matrix of single-qubit WF 

= ({as, Ag}x) = uiQx) - (s) (s; 



where we use the axis operators & = p,d,q, and A£j = — Moreover, the label 

X = S,D, linking to the two types of anticommutator defined in formulas El and 
ITT1 leads to two different covariance matrices. Nevertheless, the following results are 
independent from definition of anticommutator, and we will write the parameter X only 
when necessary. Recalling that matrix is a 3 x 3 semi-definite positive symmetric 
matrix and that diagonal elements Vu, named variances, coincide with the uncertainties 
U[£i), then the sum of diagonal elements Vu is positive, consistent withlUA 

Following the scheme of reference [32] for the continuous case, the covariance matrix 
of two qubits is built by writing formula HH1 with the enlarged set £ = (pi, d\, §i, P2, d<i, 92) 
giving a 6 x 6 matrix. A compact version of covariance matrix is given by 

A C 



V 



C T B 



(19) 



where A, B, C are 3x3 matrices. Notice that matrix elements of A and B represent the 
covariance matrix Vij relevant to qubit 1 and to qubit 2, respectively, while matrix C 
represents the inter-qubit correlations between axis operators £j and Q 2 \ Two-qubit 
covariance matrix can be easily computed by means of equations El once the WF is 
known. It is easy to show that the inter-qubit correlations Cu measure the degree of 
correlation between spin observables 0^ and . Hence their operational meaning is 
the establish the interdependence of the two constituent subsystems. At this point, the 
WF formulation of LUR's is easily achieved. Upon observing that 

u[& 91+19 g 2) ] = um + um + 2{ (g i) ® g 2) ) - (e? } ) (g 2) >} , 

the LUR relevant to the axis operators becomes 

trA + trB + 2trC > 1 , (20) 

where only diagonal elements of submatrices are involved, thus making the formula 
independent from the definition of anticommutator. This equation has the following 
interpretation: if the correlations Cu are negative and their absolute values are 
sufficiently large, than the inequality is violated and the state is entangled. This 
evidences that non-separability strongly depends on the inter-qubit correlations 
described by tr C. An important problem that deserves to be clarified is raised by those 
entangled states where correlations Cu are positive. To answer to this question, in tableEl 
we consider covariance matrices relevant to both the WF of Werner's state (including as 
well the singlet state with weight x) and the WF of the Bell state |$ + ) = (|00) + |ll))/v / 2. 
In the first case, formula I2*U1 is violated for x < 1/3, which is a correct result. Instead, 
in the second case, formula EDI is not violated, thought the Bell state is known to be 
maximally entangled. We can show that, except for the singlet state, no Bell state 
violates formula EH and the criterion does not supply information about separability. 
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The singlet case (corresponding in the table to the x — 1 case) differs from the 
other Bell states in that all the diagonal elements of matrix C are negative. The 
other Bell states have elements Cu with alternating sign, which makes the violation 
of formula [201 impossible. Nevertheless, it is clear that such correlations, thought not 
negative, contain a large amount of information on non-separability. We thus propose 
the following modified inequality as a necessary condition for separability 

z\c,A< trA+t : B -\ (2D 

i z 

whose main feature is to replace the diagonal elements of C with their absolute values. 
The effectiveness of formula |^ is confirmed by the fact that it is violated by any Bell 
states. The quantum-mechanical meaning is also clear in that, if a state is separable, 
then the absolute value of the correlations must be bounded from above. We easily prove 
that this inequality follows from equation El by resorting to the more general operator 
Q <8> / + €il <S> Q where e; = ±1. With such a choice, the LUR condition written in 
terms of covariance matrix reads trA + trB + 2 J2i > 1. To obtain inequality El it 
is sufficient to consider ti = — 1 when Cu is positive. As a final comment, we notice that 
the sum of correlations |^ thus exhibits an upper limit for separable states: entangled 
states may overcome it, and the exceeding part is an indicator of quantum correlation. 

3.5. Generalized Uncertainty Principle (GUP) and PT criterion 

In the continuous well-known separability criterion jHB E2j is obtained by 

combining the Generalized Uncertainty Principle (GUP) with the application of the 
PT criterion to the variance matrix relevant to position and momentum operators. 
For a system of two (one-dimensional) particles in a continuous space, the GUP-based 
criterion states that, if a state p is separable, one can construct a matrix M = tr(p^j) 
which is semi-definite positive under PT, namely [32] 

M = V + ^tt>0, (22) 



" J 


" 


, J = 


1 " 




_ 


J _ 




-1 



where V af3 = ({A£ a ,A£p}} is the covariance matrix, A£ a = £ a - (£ a ), £ Q = Pi} 
and [£ a ,£p] = i£l af 3 with 

n 

Following the PT criterion, given a separable state p and its WF W p , the PT generates 
a nonnegative operator p T ' 2 and a genuine WF W p r 2 still satisfying the equation l2*2l 

The extension of the previous GUP-based criterion to the discrete case requires that 
each separable state can be associated to a matrix M semi-definite positive under PT. 
Considering first the single-qubit case, we define the matrix My = \tr(p^j)] written in 
terms of the list of operators ^ = I,p,d,q. It can be easily shown that p > entails 
M > 0. The latter is equivalent to the condition 



Vf k + -e 3klX i>0, (23) 



2 
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Table 6. An example of GUP in the two-qubit case. 
M 





3 

\ 



-f o 

-? 



3 
J 


o 



M 



-t 



where covariance matrix in equation in the present case, is a 3 x 3 matrix 
and is related to the standard definition of anticommutator Condition implies 
tr(V) > 1/2, which is equivalent to the LUR equation for single qubit [2ZI- It is worth 
observing how any other choice for the set implies that M > iff p > provided 
{C,i} forms a complete basis of the space of hermitian matrices for a single-qubit. 

In the case of two qubits, once more in analogy with the continuous case, it seems 
quite natural to derive M from the set £ = (I,pi, di, q\,P2, d 2 , (fe)- Following the 
standard prescriptions [22] for calculating the GUP inequality, we have that M > 
and, equivalently, 



Cu B + i-e v (2) 

^ km ^mn ~ o^mns As 



> . (24) 

2 c mnsXs ' 

The matrix on the left-hand side is a 6 x 6 matrix that can be written in terms of 
the 3x3 matrices A, B, C appearing in equation Similar calculations show how 
M = tr(p T2 is such that M > if p T2 > 0. Then we conclude that the separability 
condition for a state p (achieved within the PT criterion when both p > and p T ' 2 > 
are satisfied) is now ensured by M > and M > 0. Notice that condition M > can 
be reduced as well to the equivalent form 



Ajk + ^tjklXl ^ Cj n 

Cu B + if y (2) 

^ km -L-'mn o c mnsAs 



>0, (25) 

2 c mnsXs ' 

where B, C are determined using once more the PT operation. Formula 1251 containing 
the axis-operator covariance matrix is the core of the two-qubit PT criterion. In table 
IH1 we illustrate the application of the present criterion to the Werner state. In this case 
V = M and the eigenvalues of matrix M (relevant to p T2 ) are positive for x > 1 (rather 
than for x > 1/3). Unfortunately, this means that the GUP is not violated so that 
the criterion does not give information about separability. This can be explained with 
the fact that, when using the set of operators (pi, di, qi,P2, d 2 , 52), the nonnegativity 
of M is only a necessary condition for separability. In order to cure this problem we 
have generalized the GUP-based criterion by using the enlarged set of operators £j <8> 
which leads to a 9 x 9 matrix M. In this case, we could have a violation of the positivity 
condition under partial transposition. This result will be discussed in a separate paper. 
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4. Conclusions 

In the present work, we have considered the WF defined in pQ focusing our attention on 
two properties of the two-qubit WF, the negativity and the covariance matrix, which 
are useful in the characterization of entanglement. After reformulating/generalizing the 
PT, LUR, and GUP-based separability criteria in the WF formalism, we have tried to 
evidence what features of the WF and of its covariance matrix are able to reveal the 
presence of entanglement. 

In section I^TTI we have found that a two-qubit WF relevant to a separable state can 
not assume values lower than (1 — \/8)/4. In section EOl we have recast the PT criterion 
ITU in terms of WF by means of inner-product rule^J Based on this result, in section EHU 
we have shown that the non-separability of p entails the presence of negative elements 
in W p or in W„t 2 . Interestingly, these facts relate the main non-classical feature of the 
WF (the presence of negative values) to the presence of entanglement in the two-qubit 
system. Considering the separability problem within the LUR criterion, in section EPl we 
have reformulated it in terms of covariance-matrix elements of WF|2I1 in particular, we 
have found a stronger version of the LUR criterion (illustrated by formula l2l"j) once more 
involving the covariance matrix. This generalized criterion, which has been tested both 
on Bell states and on Werner states, evidences that the presence of strong correlations 
can be used to detect non-separability. Finally, in section IH.5[ we have studied the 
analogue of the GUP-based separability criterion (continuous case) from the viewpoint 
of discrete WF's. We have shown that adopting the same procedure of the continuous 
case leads to criterion [213 The latter does not succeed in detecting entanglement as 
a consequence of the fact that the set of operators used to build the discrete GUP is 
too small. In order to cure this problem, we have enlarged such an operator set thus 
obtaining that M > p Ti > 0. Such an equivalence provides the basis to extend 
in an effective way the GUP-based separability criterion from the continuous to the 
discrete case. 

Future work about entanglement properties of the two-qubit WF will be developed 
in two directions. Our first objective is to derive, relying on equation [213 the explicit 
form of a generalized GUP-based separability criterion from a suitably enlarged operator 
set. A second important problem which deserves to be deepen is to establish how the 
presence of negative elements in WF's W p (and W p T 2 ) relevant to entangled states is 
related to the violation of inequality E] issued from the LUR condition. Such aspects 
will be investigated in a separate paper. 
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